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- technique HEATEST (HEATing ESTimation) are made. The 
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with changes in time step, thermocouple length, and number 


of discrete node points. The estimation of the convective 
heat transfer coefficient and thermal conductivity were 
excellent, 
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HEATING PARAMETER ESTIMATION USING COAXIAL 
THERMOCOUPLE GAUGES IN WIND TUNNEL 


TEST ARTICLES 


I. INTRODUCTION 


1.1 Background 

The determination of heat transfer rates on hypersonic 
configurations upon reentry is important for the survival of 
the vehicle. The problem is that the heat rate is not a 
quantity which may be directly scaled from model tests in 
wind tunnels. However, the parameters which make up the 
heat rate equation (thermal conductivity and specific heat, 
for example) can be scaled from which the heat rate may then 
be calculated. Wind tunnel heat transfer measurements have 
traditionally used a thin walled test model fabricated with 
thermocouples mounted on the inside skin surface. The 
"thermal model" then yields a heat rate based on temperature 
measurement from the thermocouple. Another technique uses a 
coaxial thermocouple gage mounted in a thick skin model, A 


discussion of the two methods follows. 


The Traditional Thin Skin Model 


The "thermal model" of a traditional wind tunnel thin 


skin model assumes that all of the heat penetrates the thin 
skin via conduction to a standard thermocouple gage mounted 
on the back face. No lateral conduction is assumed and 
since the emittance of the steel model is low, radiation is 
assumed negligible. A typical heat rate measurement data 
point is aquired by injecting the cooled model into the wind 
tunnel at a known temperature and time and by measuring the 
temperature at later times. The model is then removed from 
the tunnel, cooled, anda change in configuration is made 
in preparation for the next injection and subsequent dat 

point. There are three very severe limitations associater 
with this technique (Ref 1).The first is the inability to 
acquire more than one data point during any one injection. 
The thermal model simply does not allow for the type of 
change in configuration or model attitude which can be 
accomplished using dynamic testing techniques (to be 
discussed later). Associated with this limitation is the 
long cooling time between each test which significantly 
increases the cost per data point for the overall test. A 
second limitation is the special thin skin model which must 
be fabricated, further contributing to increased test cost. 
Finally, the assumption of no lateral conduction through the 
model may in fact be a poor assumption at some critical 
locations with large curvature. An alternate type of gage, 


the coaxial thermocouple, can eliminate these limitations 


with an overall effect of reducing time and cost. 


\ 


A New Application For An Old Thermocouple 

The coaxial thermocouple gage is shown in Fig. 1.1. It 
consists of a constantan (a metal alloy) jacket surrounding 
a chromel core with a thin layer of insulation separating 
the two metals. The coaxial gage is mounted ina steel 
model thick enough so that the thermal pulse is not sensed 
on the backface (ie. the model wall is considered a semi- 
infinite slab). The thermocouple surface is formed when the 
gage is lightly sanded to match the contour of the model. 
Some gages are avaibabie with backface temperature 
monitoring to assure that the thermal pulse does not reach 
the backface in any given run so that an analytical 
integration of the heat equatior can be used to determine 
the heating rate history, The backface temperature 
information prior to this investigation is not factored into 
the data reduction process, however. 

Operation of the coax gage is based on uniform 
conduction along the gage length which would necessitate the 
model be made of a material with similar thermal properties 
(Ref 3). The thermal properties of stainless steel match 
very closely with the gage properties, therefore, the 
presence of the gage is negligible. The matching of thermal 
properties also enhances accuracy. A cOaxial gage which is 
matched thermally with the model allows an isothermal wall 
assumption, whereas other gauges such as calorimeters and 
thin film gages are not thermally matched, and cause a non- 


isothermal wall. Measured heat transfer can be in error by 


way Ye 


{0} 


(S,) 


{A]{U} + {b} + W(t) 


[A](Sy)} + (dy) (3-1) 


These equations are solved using a tridiagonal algorithm in 
subroutine TPS3 for the temperature states, and subroutine 
SENS for the sensitivity of the temperature to the kth 
parameter. Propogation of the covariance, P, of the 
temperature state at each node is accomplished by the 


approximate difference equation, 


P(t.~) = 6(At) P(t, 47) st (at) 
th 


+ 6(t -2)Q6T(t -A) dA 


si (3-2) 


where ¢ is the transition matrix and where the - and + 
superscripts are used to denote the expected values before 
an update (or a priori) and updated (or a posteriori) 
values, respectively. This calculation is made in 
subroutine TPSOSP2. 

A model of the temperature measurement process must be 
used for the Kalman filter equations. The measurement 


equation to identify thermocouple location is, 

Y(t,) = H (U(t,)} + {pg} (3-3) 
The updated temperature is calculated by, 

B(t,*) 


n ) = U(t, ) + GE(t,) (3-4) 
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time. They are found by employing a Kalman filter - a set 
of recursive equations that optimally combine the 
propogation of the model equations with measurement updates 
at each sample time. After the entire temperature - time 
(state) history has been calculated, a gradient algorithm is 
used to solve for best estimates of the parameters according 
to a maximum likelihood criterion. 

The second type of estimate consists of the parameters 
defined in the parameter vector,{@}, as given in Equation 
2-10. These parameters remain essentially constant 
throughout a transient maneuver profile such as a pitch 
sweep and are estimated based on data from the entire 
maneuver history. The process of estimating states and 
parameters is then iterated for convergence to some optimal 
estimate. 

The method used to estimate the states and parameters 
is formulated from stochastic estimation theory and is known 
as adaptive estimation. A detailed development of the 
estimation equations is beyond the scope of this thesis and 
the reader is referred to References 7 and 8 for more 
detail. 

The thermal model equations for temperature and 
sensitivity have already been written in the matrix 


stochastic estimation form of Equation 2-12 as, 
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Figure 3.1 HEATEST Algorithm Summary 
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III. HEATEST OVERVIEW 


The HEATEST program was originally developed to deter- 
mine and model heat rates from the Space Shuttle Orbiter 
thermocouple data, hence most of its nomenclature references 
flight data samples and trajectory samples. The wind tunnel 
equivalence of the trajectory sample would be the test 
section conditions at the time of the sample (ie. density, 
velocity, pressure, etc.). The flight data are the thermo- 
couple measurements from the coaxial gages. 

An algorithm summary of the HEATEST program is given in 
Figure 3.1. The initial conditions for the temperature 
profile, U(t), and the initial covariance, P(t), are 
specified at the start of the wind tunnel test run. Heating 
model initial parameters, and the initial reference values 
for the heating model are read in as inputs to the program. 
Initial sensitivities of the state are specified to be zero. 
The node point structure throughout the depth of the 
thermocouple is then calculated from the input of the length 
of the gage and the number of node points. 

Two types of estimates must be made in order to 
describe the thermodynamic environment in the wind tunnel, 
The first type are the state estimates, which are defined by 
each node temperature. These state estimates are not 
constant since the temperature varies throughout the 


maneuver, and hence, must be estimated at each node point in 
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The {b} vector is not used directly for the covariance 


equation, but is approximated by an error model given by, 


Q (Bhardreg (Tay, ) Ax/6,k]? (2-15) 


error ~ 
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ee where S is the sensitivity. The i subscript identifies the 
7 mode point and the second subscript identifies the 
particular parameter number. The sensitivity equations may 


also be written in the familiar forn, 
[A'](S, 47} + {4} = 0 (2-12) 


The sensitivity equations are developed and summarized in 


Appendix B. 


Propogation of the covariance of the temperature state 


q at each node requires the equations to be of the form 
¥ (U} = [A](U) + {b) 
ae and (S; 4} = [A]{S,} + Cay) (2-13) 


Substituting the definitions of Equations (A~14) into 
Equations (2-7) and (2-8) and rearranging yields a common 
tridiagonal [A] matrix for the above equations which is 


shown presently, 


Gia ead COLLAB bal ited.) Fe) : 
1 1 


[A] = RM -(RM,+RP, ) RP 
get, : ret, 


(2-14) 


(A']{0,;") + (b} = 0 


where [A] is an n xo tridiagonal matrix of material 


properties and {U,") is the n-dimensional column vector of 


unknown temperature at each node point for each time. 


In general, the unknown parameters in this model 
formulation are the heat transfer coefficient intercept, hg, 


the slopes hoi and hao > and the scaling parameters for 


specific heat and thermal conductivity, $4 and by respec-— 


c 


tively. 


These parameters may be defined as a vector, , of 
unknown parameters for use in the system identification 


scheme as, 


d= (ho, bays bags be Oy} 7 (2-10) : 


The primary purpose of the heating estimation program is to 


obtain best estimates of these parameters during transient 


test maneuvers. To estimate these parameters it is 
necessary to calculate the model sensitivity to each unknown 


parameter. 


2.2 Sensitivity Equations 
The derivative of Equation (2-7) with respect to each 
parameter yields equations of the same form as Equation (2- 


9) from which the HEATEST program propogates the 


For the 


sensitivity. example, sensitivity of the 


temperature with respect to ho would be written as follows, 


aU = aU = Sho = Si1 
(2-11) 


level defined by the superscript n, 


(Ua etsy 4 Za (U,2»%)4 + 4(0,2> *)3(0,2-s*1 = U,7°*] 


-3(0,7"5)4 + 4(u 829570, eet (2-6) 


Substituting Equation (2-6) into (2-5) yields, 


Ny a-1 = 
p6,cAx Ui U, 6yk141/2 ‘ $4%141/2 n 
ee ee ee ee — a ae ee eI | | a ere 
2 At Ax Ax 


- eo f-3(U, 7 §)444(0, 2 $) 30,2 842-048) 4] 
+[ both. y(a-a,)+hyo(a-a5) I Begg (Tg y—Uy 8?) 


(2-7) 


The model equation for the interior node points, 


(i=2,imax), yields, 


pél icAx Fig alae 6yki-4/2 


sacar | 

At Ax 

6, (Ki-1/2*Ki41/2) 
oo U; 

Ax 

64k i+3/2 

Lene os | 
Ax 


(2-8) 


Equations (2-7) and (2-8) can be rearranged into the 


familiar matrix forn, 


‘, s, 4, 4, 4- 


Poy 
sits 


Ps 2 
oe La <i igen 


parameters other than those included in the reference heat 
transfer coefficient are summarized by the static transfer 


relation or heat transfer coefficient ratio, h =h/h 


bar ref: 


Here, the ratio is assumed to be piecewise linear with 
respect to angle of attack as derived from Lagrange 


Interpolation Theory (Ref 6 and Ref therin). 


boar = [hg + bgy(a - ay) + bggla - ag)] (2-4) 


where hg is the magnitude of the heat transfer coefficient, 
h at the reference conditions ,a, , specified by the one 
subscript. The hesting parameters ho » boy - and ha? are 
considered to be unknown and constant over a prescribed time 
period, and will be estimated by the HEATEST program. The 
parameters correspond to a derivative with respect to 
deflection angle of the model. Thus, for constant step 


size, the model equation at the surface node, (i=1), 


becomes, 
n_ n-1 = 
a les A a cl Sa Se UE 
2 At Ax 2 


6yk141/2 
+ 


U9 ~ eof(Uy*)* - (U42) 4] 
Ax 


+ [hothgy (a-ay) thay (a-an) Ibe eg (Tgy- 04") (2-5) 


The non-linear radiation term is quasi-linearized on an 


iteration level defined by the superscript s and by the time 


e radiative emissivity 

o StefanBoltzmann constant 

c material Specific Heat 

p material density 

k Thermal Conductivity 

6, Specific Heat scaling parameter 

6, Thermal Conductivity scaling parameter 


where 


The material specific heat and thermal conductivity are both 


scaled by the two factors 6, and by: respectively, hence the 
value for c and k will remain unchanged. The parameters 6. 
and ¢, will be estimated by the HEATEST program. 
Coefficients with subscripts which are less than one or 
greater than n are zero. The radiation and heat rate terms 
are also zero except at the surface node, Equation 2-1 
includes terms due to conduction from adjacent node points 
Ki-y/2/4xj-) » surface radiation eot;4 , and the convective 
transfer of energy as obtained from the heating model. The 
resulting system of implicit difference equations must be 
solved simultaneously. 

The heating model for the convective transfer of energy 


is based on Newton's Law of Cooling, 
i (Ty - T) (2-2) 


Non-dimensionalizing by a reference heat transfer 


coefficient, hoes yields, 


4 = Byardreg (Tay - Uy) (2-3) 


The dependance of the heat transfer coefficient on 


2.1 Temperature Equations 
Across section of the one-dimensional model is given 


in Figure 1.1 and below as a typical coaxial thermocouple 


gage. 
Grad oP 
Foi ft fe 
q = 0 
6 dx dx 
i-1 i-1/2 i i+1/2 i+1 
Geonv 


An energy balance is performed on each element. The thermal 
conductivity, k, is taken as an average between each node. 
Fourier’s Law of Heat Conduction throughout the gage, the 
Stefan-Boltzmann Law for radiation and Newton's Law of 
Cooling for convection on the surface face yield a system of 


mn nonlinear differential equations of the form: 
(ps8 .cgAxgtp 346 ,¢5-yAxy-y)/21((0,9-U,2 1 )at] 
= Ohya 2/hxq-4)054" 
Oy Uk y-4/2/Axg-athings2/Axiy 1104” 
+ Oy (kia gs2/ Arie IU 544” 


oa 4_ 4 Zs ze, 
eo(U; Ui-y thy arbres (Tay U;) (2~1) 


@age using a Kalman filter(Ref 2). The purpose of this 
investigation is to incorporate the appropriate thermal 
model equations into the HEATEST program for application to 
coaxial thermocouple gages as used in the wind tunnel. Two 
reasons for doing this are, 1) replace the analytical with 
the semi-infinite assumption, ie. extend the run time or 
shorten the gage, and, 2) estimate thermal properties. 
Testing and verification of the modified program is 
necessary for verification of the validity of the results. 
Simulated data are generated by an analytical solution, and 
are processed for testing purposes for which the results are 


known. 


1.3 Overview 

A development of the pertinent temperature, sensitivi- 
ty, and covariance equations will be developed for introduc- 
tion into the HEATEST program in Chapter II. Details in 
format for programming may be found in Appendix A & B. 
Chapter III is an overview of the HEATEST algorithm and 
shows how the equations developed in Chapter II are 
utilized. Chapter IV outlines the method for testing the 
program and offers a discussion of the test cases made and 
results. Finally, conclusions about the validity of the 


modifications, and suggestions for further improvement are 


included in Chapter VY. 


up to 40% because of the non-isothermal (Ref 2). The same 
rugged model built for pressure measurements can be used for 
temperature measurements as well, which would further reduce 
the wind tunnel costs. The data reduction technique, which 
uses the temperature time history, eliminates the 
requirement to fix the model configuration or attitude 
during any one run, hence dynamic testing techniques may be 
used similar to the flight test technique used for the Space 
Shuttle Orbiter (Ref 4). The model may be swept in angle of 
attack, for example, to determine heat rates as a function 
of angle of attack. All of these attributes along with a 
short response time and no required calibration (ref Knox) 
combine to yield the wind tunnel engineer a tool of marked 


improvement over previous methods. 


1.2 OBJECTIVES 

A method of analysis to identify the aerothermodynamic 
flight environment and update the thermal model of the 
ergineering simulation of the Space Shuttle Orbiter was 
designed by the Air Force Flight Test Center. This method 
is in the form of a digital computer program called HEATEST 
(HEATing ESTimation). The program provides a correlation of 
the heating as well as a heat rate time history. The program 
integrates numerically instead of relying on some analytical 
assumption. It satisfies a maximum likelihood criteria for 


each parameter and obtains best estimates for the 


temperature at discrete nodes throughout the length of the 
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where 


Q 
WW 


P(t, )HT(BP(t ~1)H7 + RLJ~2 


w 
i 


Y(t,) - HU(t, ) 
The updated sensitivities are calculated by, 
S,(t,*) = [I-GH]P(t,)(I-GHI]T + or,GT (3-5) 
The updated covariance is calculated by, 
+) = [I-GHIP(t, )[I-GH]? + oRr,G? (3-6) 


To alleviate the problem of imprecise initial 
conditions, a fixed point smoothing algorithm has been added 
to the HEATEST program. Details of the smoother and its 
effects in the adaptive estimation scheme may be found in 
Reference 6. 

Finally, the best estimates of the parameters are then 
estimated at the end of a specified time segment by the 


gradient algorithm, 


yom ae {a7F/a02]-lar/ae = @ + s71s (3-7) 


where, 
. ~)qT 
ij =D Sie s )HT(HP(t,~)HT+RmJHS, y(t.) 
n= 


N 
Sy So Sipe ig Je TEMS TUG Snes 1) 
n= 


The matrix J is an approximation for the Jacobian or 


conditional information matrix and is given in component 


form by Jijzn + The score vector, Sy, is used to approxi- 


mate the gradient of the likelihood function for a large 
mumber of time samples. 

Using these equations, best estimates for the 
temperature time history (states) at each node can be found, 
as well as the deviation in temperature as provided by the 
covariance matrix. Also, an estimate of the parameter 
uncertainty is provided by the Cramer-Rao bound. The 
Cramer-Rao bound relates the conditional information matrix 


to the covariance of the parameter estimate. 
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To test the validity of the program modifications, a 


set of contrived data was generated. It’s development 
assumes that the heat rate due to convection at the surface 
node is constant and equal to the heat rate due to 
conduction at the surface. The heat rate due to convection 


is given by, 


or 


2 = Boarkref(Tay ~ Ty) (4-1) 


where hha, is defined as in Equation A-12. The equation for 
the heat rate due to conduction assuming a one-dimensional, 


homogenous, semi-infinite solid is as follows (Ref 9), 


t 
(pek)2/2 ( at(c) = at 


q = so 
n dt (t-7) 1/2 
re) 
(4-2) 
where t = timefromstartofheating 
T(t) surface temperature rise 


dummy variable of integration 


a 
tow 


Equating Equations 4-1 and 4-2 yields, 
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Sa 


(pek)1/2 [a dt 


a= bpgrtser(Tay Ty) * _—__— 
ar’re aw ‘Ww : ae feng il? 


(4-3) 


Two different expressions for the derivative of the 
wall temperature with respect to time were used. The first 
implied a linear change in temperature with respect to time 
yielding a constant for dT /dt and the second expression 


assumes that temperature was a quadratic function of time as 


shown, 
Linear Quadratic 
Ty, = bt + T, = at? + bt +c 
oe = bd ae = 2at + b 


Solving Equation 4-3 for T,, so that h and q are constant 


and after making the indicated substitutions and 
integrating yields, 


Doarbrer 
(4-4) 
Quadratic assumption 
2 pekt fb + 4at 
Taw = Ty + ) 
Dharbrer n 3 
(4-5) 
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A short computer program was written to produce a 
temperature-time history in the data tape format for the 
HEATEST progran. For the above equations, hief and hbpa, 
“ (the estimated parameter) were set equal to 1 and the 
coefficients a, b, and c, were selected to yield reasonable 


values for Ty: 


The reference test case was taken to be a 10 sec. 
simulated wind tunnel test run using the linear data pro- 
vided from the previous section. Thermocouple samples and 
test section samples were provided at the rate of one sample 
per second. The objective was to examine the rate of 
convergence of the temperature states and to estimate the 
te first parameter, hyp. Recall from Section 4.1 that the data 

was generated to yield a value of one for hg- Also of 
interest, was the validity of the model to the semi-infinite 
solid assumption (ie. no change in the temperature at the 
back face node throughout any specified time segment). 

The input data is shown in Figure 4.1 and is generated 
digitally depending upon a desired time step (del t). The 
initial temperature throughout each gage is assigned a value 
of 60°F. Figure 4.2 shows the input temperature values for 
a del t = 1 sec. as used in the reference test case, 

Figure 4.3 identifies the temperature state at the 2 
sec.(lowest curve), 6 sec.(middle curve), and 10 sec.(top 


curve) times following the first iteration through the 
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updated HEATEST progran. It clearly shows that at the back 


face node, the semi-infinite solid assumption used to derive 
the data is violated. This is indicated by the change in 
backface (node 6) temperature with time. Note also, 
however, that the temperature gradient at the back nodes 
(between nodes 5 and 6) is zero due to the adiabatic wall 
assumption. It should also be pointed out that surface node 
temperature response to the given input was immediate with 
no time lag. 

Figure 4.4 is similar to Figure 4.3 except the 
temperature states and parameter estimates have been 
iterated to convergence, in this case, three times. 
Overlaying the two figures shows no perceptible difference 
between them, and the data shows no variations in values 
until after the second decimal point. In spite of the 
response of the back face node which would ordinarily 
invalidate the test, the estimated value for hg was .99598, 
within .4% of the desired value of 1! Several test cases 
will be compared to this reference by examining changes in 
time step, thermocouple length, and the number of node 
points. Also, an examination of the ability of the program 


to estimate the other parameters, follows. 


4.2.1 Changes in Time Step 
Figures 4.4, 4.5, and 4.6 show temperature states at 2, 
6, and 10 sec. for del t equal to 1, .5, and .25 sec., 


respectively. All three curves required three iterations 
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— 


eS for convergence. The curves are all similar in shape with 


almost no perceptible differences. However, if the 10 sec. 
curves from del t = 1 sec. and del t = .25 sec. are 
overlayed as in Figure 4.7, a small difference may be noted 
at the backface nodes indicating that, indeed, a decrease in 
time step will yield a profile which will more closely 


approximate the model. 


4.2.2 Changes in Thermocouple Length 

Changes in thermocouple length offered the most 
dramatic changes in temperature state as can be seen in 
Figures 4.8, 4.9, and 4.10 where the lengths range from .1 
ft., .OSft., and .025 £ft., respectively. The two longer 
lengths converged within three iterations while the short 
thermocouple length took four iterations to converge. The 
extra iteration is most likely due to the large deviation 
from the model and the large differences in temperature 
state from one time step to the next. Figure 4.11 compares 
the temperature state of each thermocouple length after the 
10 sec. run. The lowest curve is associated with the 
longest thermocouple, and the upper curve is associated with 


the short gage. 


4.2.3 Changes in Number of Node Points 
For this comparison, a 10 sec. run with 6 node points 
and s 3 sec. run with 6 node points (Figures 4.12 and 4.13) 


run with 12 node points and 


will be compared with a 10 sec. 


a3 sec, run with 12 node points (Figures 4.14 and 4.15). 


Each of the test runs were converged by the third iteration. 
At both of the different run times, increasing the number of 
nodes yielded a solution which more closely approximates the 
model (ie. the semi-infinite solid at the back face node) 
and gave correspondingly better estimates for ho- The 
comparison of temperature states is better represented by 
Figures 4.16 and 4.17 which directly compares 6 nodes and 12 
nodes interspersed evenly throughout the .05 ft. long 


thermocouple for 10 sec. and 3 sec. run times, respectively. 


4.2.4 Parameter Estimation 

The ability of the algorithm to estimate hy has already 
been discussed. To summarize, even when the output 
temperature states clearly violate the semi-infinite solid 
assumption used in generating the data, the estimated values 
of hy remain within 15%. The 15% is a worst case number 
derived from the short thermocouple using a course grid for 
a long run time. An average deviation which considers all 
of the test cases evaluated is closer to 3%. 

To estimate é,, an erroneous data value was given for 
ee Kgata’ Wherin the program iterated to a value for $6, which, 
when multiplied by the erroneous Kaata’ Would yield the 


*. correct value, ie. 6, Kgata = Kcorrect ° The erroneous 


Kaata Which was input into the program was 14% in error of 


the Koorrect Valine: The value provided by the program for 


$, when multiplied by Kyat, yielded a value within .9% of 


vy 


value! 


the Ecorrect 
The estimate of 6, was not as successful, however, 
After 6 iterations, the solution was diverging from the 
expected value. A suspected sign error in the 6, 
sensitivity calculation is the most likely cause. 

The quadratic input data was used as a comparison to 
the linear data to challenge the algorithm, ic. the more 
complicated the input the more difficult the estimation 
process. No direct comparison may be made of temperature, 
however, since the input data is different. The parameter 
estimation of hg using the quadratic input data was still 


excellent yielding .1%, while the estimate using linear data 


was somewhat better at .05%. 
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the utility of the program to be used in wind tunnel runs of 


much longer duration. 


5.2 Recommendations 

Actual wind tunnel test data from coaxial gages needs 
to be analyzed by the program to instill more confidence in 
the results. This would require that the data tape format 
be modified to be compatible with the program inputs. 

Also, a prescribed model for Loar as a function of 
angle of attack needs to be input to determine the ability 
of the program to estimate the piecewise linear derivitives, 
hed and h,o. The additional input data would then be the 
wind tunnel model angle of attack at each thermocouple 
sample time. 

Another potential modification to the program would be 
to use a second order time derivitive approximation as 
opposed to the current first order approximation. It is 
suspected that the increased accuracy would improve the 
state estimates particularly for large time gaps in 
thermocouple data. 

The temperature state estimates might be improved by 
incorporating a variable grid. An exponential grid 
generation scheme would concentrate node points near the 
surface where the largest temperature gradients exist. 

A final recommendation, would be to make the program 
more ‘user friendly’ and to publish documentation much like 


aousers manual. 
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¥, CONCLUSIONS AND RECOMMENDATIONS 


5.1 Conclusions 

The program works for this simplified test case of 
known parameters and number of node points. The ability of 
the program to determine the temperature state appeared very 
good but definite conclusions cannot be drawn without 
comparing this data to actual thermcouple data. The ability 
of the program to estimate the parameters, hg and ¢, was 
demonstrated extremely well. Re-examination of the 6, 
equations must be made before further conclusions may be 
made about the ability of the program to estimate this 
parameter. It is theoretically possible, however. 

The objective, to validate the use of the modified 
HEATEST program for use of coaxial thermocouple gages on 
wind tunnel test articles, has been met for the special case 
of this analytical model. The program is capable of 
determining temperature states and estimating parameters 
with a high degree of accuracy. A note concerning the semi- 
infinite slab assumption needs to be emphasized. This 
assumption was made only for the data generation program to 
yield data for which the analytical solution was known. As 
mentioned in Chapter 1, some coaxial gages are available 
with backface temperature monitoring which would also 


provide another temperature measurement to enhance 


estimation of 6, and 6,. This feature would greatly enhance 
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APPENDIX A 
‘ Derivation of equations using a one-dimensional energy 
balance formulation are given as follows, 
‘ 
Grad Ss 
as 
: q = 0 
; a dx dx 
4 i-1l i-1/2 i i+1/2 i+1 
Geonv 
r Energy in the left face = - kdT/dx = a, 
Energy generated within the element = q dx = 0 
Change in internal energy = pe(dT/dar)dx 
"7 Energy out of right face = ~ k(9T/8x) 44, 
= -[kaT + 0 skaT\dx] 
ax ax ss) 
r Then, combining the above and using Fourier's Law of Heat 
< Conduction, ie, 
energy in , energy within = change in + energy out 
e left face the element internal energy right face 
yields, 
- koT 4, qdx = pcdTdx -fkaT , a ‘ex a 
e ax at ax O@x\ dx (A-1) 
- or, 
peat _ ieee 
e e at Ox\ ax (A-2) 
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<i 
ase 


or, replacing T by U, 


pcu, = (kU,), (A-3) 
or, 
k, k: 
i-1/2 i+1/2 
pcU,y = a (U,;_,7-0;") TS (U;"-U;,,7) — 
Axi-4/2 Axi41/2 Ax 
(A-4) 


where AX may be written as, 


Ax = AXiaas2 + Axj-1/2 
2 


then, writing the time gradient in first order backward 


difference form and expanding yields, 


na. 
Uye-Uz-4" = 2Ki-1/2 
a Se eg 
a Axi-a/2(Axjea/2t4x5-1/2) 


pe 


k;_ 
7 2ki-1/2 u,;* 


AX 5@4/2(AX441/2t4%4-1/2) 


2ki+3/2 
s —_—________________ 0,” 
AX i44/2(A%541/2t4%4-1/2) 


2ki+1/2 
+ ——————_—_______________ 0, ,” 
Axign/2(Axj4a/2tAxj-1/2) (A-5) 


then, specifying equal spacing for each node point, 
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-!. 


oe 


AXigas2 = Axj-1/2 


and the equation becomes, 


nw n-1l . 
pe 


Ki-as2 
At Ax? 


no 
Uj-4 


k, k; k: 

~1/2 +1/2 +1/2 

rm i-1/ ae. 1/ vj. + it+1/ U n 
Ax? Ax? 


(A-6) 


Now, instead of estimating c and k directly, define two 
scaling parameters $6, and ¢, such that c and k will remain 
constant. These two parameters are estimated by the HEATEST 


program. 


n_ ~1 
Vie Ua ki-1/2 : 
pi,cAx ———————-_ = 6, ————— 0-1 

At Ax 

k ; + & k; 
i-1/2 i+1/2 i+1/2 n 

- 6y(———- s,s ———— ) 0," + 6, iggy 

Ax Ax Ax 


(A-7) 


This equation is applicable at all interior (i # 1,i # Loex? 
points. 

For the back face, assuming a semi-infinite solid, i = 
and the equation becomes, 


Imax 


n_y n-1 
2 At Ax Ax 


(A-8) 
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a 2 For the front face, (i = 1), the effects of radiation away 
from and convection toward the solid surface must be 


accounted for. 


The radiation is modeled using the Stefan-Boltzmann 


Law, 

& cs 4 _ 4 - 
on q = eo(U, Up ) (A-9) 
a where g radiative emissivity 

Ps o Stefan-Boltzmann constant 

: U Temperature (°R) 


The convective transfer of energy is modeled using 


Newton's Law of Cooling, 


a 7 


qa = K(T,, - Ty) (A~10) 


Non-dimensionalizing by a reference heat transfer 


coefficient, h,.¢ yields, 


9 = Byardree (Tay ~ Uy) (A-11) 


where, hpar™ convective heat transfer coefficient ratio 


av = adiabatic wall temp of test article 


The dependance of the heat transfer coefficient on 
parameters other than those included in the reference heat 
transfer coefficient are summarized by the static transfer 


relation or heat transfer coefficient ratio, h/h. og. Here, 


the ratio is assumed to be piecewise linear with respect to 
ve. angle of attack as derived from Lagrange Interpolation 
e . Theory (Ref 6), 
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~ tool, : 
i Pir aeene 


= 4 


Byar = D/byeg = (hothgy(a-a4)t+hgo(a-ay)] (A-12) 


where hg is the magnitude of the heat transfer coefficient, 
h, at the reference condition specified at a,. Combining 
Equations A-7, A-9, A-11, and A-12 and evaluating at node 
one yields, 

poochs Uy ~ Uy™ 2 = —bykgsi/2 
2 At Ax 


+ 6,k U . 
ee got 084 = (UEDA 
Ax 


+ Cho t+hgy (0-04) thay (a-ag) hyeg (Tay-U3”) 
(A-13) 


Using the quasi-~linearization as developed in Equation 2-6, 
the resultant form for determining the temperature time 
history at each node point is given in Equations 2-7 and 2- 


8. 


The matrix form for the equations may be found after 


defining the following, 


RCX, RCX; 
RCX; = pé.cAx RCX, = —— RCX; = 
2 Zz 
6yki-1/2 
Ax 
Oxki+i/2 
RP RP; = 0 
Ax 
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Fe aay ae ae ee Oe a Ce ee en tite, We Crees 


ct 
{ oe soe RCX, 
a Jae BBB, = a RM, +RP, + 4¢0(0,7"*)3 + bya bres 
as 
ne RCX, 
. BBB; = + RM; + RP; 
o At 
= (A-14) 
' 
- 
o Then, using Equations A-14 in Equations 2-7 and 2-8 yields 
= the matrix form of Equation 2-9, 
. 
[A]{U;"} + {b} = 0 (2-9) 
a where, 
‘* -1 RP, / BBB; 0 
RM, / BBB, 
3 [A] = 
‘ RP; / BBB; 
0 RM;/BBB,; -1 
a (A-15) 
Ps and, 
Bai n,8s)4 4 n-J 
: eo[(3U, yt Des D+By arbre g Tay tRCXy0y 
At 
2 BBB, 
(o} =) Rex,u,"? 
ee At / BBB, 
|e 
Rcx;uU;" ? 
At / BBB, 
(A-16) 
e ® 
= $1 
@ 
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APPENDIX B 


The derivation of the sensitivity equations. The 
derivative of Equation 2-7 with respect to each parameter 
yields equations from which the HEATEST program propogates 
the sensitivity. A vector of parameters is formed and the 


sensitivity notation is as shown, 


# T = 
e _ {ho, hot hg2> bo 6,)) Sik = U 


k 


L—— parameter no 
Lio node point 


Defining, bya, = [ho+hgy(a-a1)+(hgo(a-a2)] , the sensitivity 


equations at node one are, 


n_ =1 = 
pA chx $1,4"~8; 4” 64 (Ki4i/2) 


e)5— ——— S Se ed 
2 At Ax 


n 


by (Ki+1/2) 

= ny)3 n _y. 2 

a a Soi4° 4e0(U,") Si, + Deeg (Taw U1 ) 
x 


n 
$151 href bar 


pBocAx Sy 27-Sy 272 by (Ei4n/2) 


9%: {= $12 
2 At Ax 


by (ki41/2) 


aa 3 n 
S. 2” 4eo(U,") $12 
Ax 


ma - n 
+ (a al)hy og (Tyy-U,") $142 heefhpar 
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PO cAx Sy 37-Sy 52 2-64 (iggy) n 
6s6 = ——— $13 
2 At Ax 
6, (ki+1/2) 
2 a a 4e0(0,")3 $13” 
Ax 
+ (a-a2) hy o¢(Tyy-U17) az $1 3 bret hbar 
(B-3) 
-1 
ee cd 4°7 81,4" ~2PRKi+4/2 ge te SPE ek 
i = ——~— 81 4 
At pb cAx pb. 7cax? 
26+k 26,k 
: byki41/2 5 date 6yk141/2 2 
p$ochx? ¢ pb *cAx? 
8eo(U,")3 2eo(U, 4-044) 
Sa Bik Pe ee 
pb cAx p6,2cAx : 
2bpartref 5, ,2 2hyarbref(TayV) 
Ae ge et hip 
pé.cAx pb, 7cAx 
(B-4) 
pbocAx Sy 5™-Sy p72 by kaa 2 a 
e5:  ——— 2. ne BS 
2 At Ax 
k k 
_ StH/2 ony beki41/2 5. 8 4 .2ti/2 4 a 
a Tyee 2,5 2 
kx Ax Ax 


i 3 a 
4eoU,"S1 5 - Sy shrethpar 
(B-5) 


The sensitivity equations at the interior node points are as 


$3 


7 . * 
aera ee Sie Tat 


2 eS ew 


follows, 


®; , i= 1, 2, 3 


Sin Sik 6yki-1/2 
p6.cAx CO FO Si-a,k 
At Ax 


6% (Ki-4/2+K%i41/2) 


i,k 
Ax 
$y Ki44/2 
‘i it+1,k 
Ax 
(B-6) 
@4 
Sing Si,g™ } byki-as2 nn ?k¥i-1/2 z 
Cee ge Rede, eee ad 
At p6.cAx pb, “cAx 
$y (ki-1/2*Ei+1/2) b 
Fi erin a 
p$ocAx Ax 
fy (ki-a/2*Bieas2) | Oe Ki+a/2 
eG aaa Ne 
2 2 » 
p,“cAx Ax ph cAx 
6yki+1/2 * 
Se Cia 
pb 2cAx? : 
(B-7) 
Os: 
2 -1 
Si,5° Si,5" byki-4/2 
pO.cAx —_———— nt S47 Og 
At Ax 
$e (kK i-1/2+Ki41/2) 
Z aes 


Ax 
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6xki+s/2 Kki-1/2 
for Ses 
Ax Ax 


(Ki-a/2 + ki4+1/2) é ki+1/2 
. + a 


Ax Ax 


(B-8) 


The backface equations are of the same form as the node 1 
equations without the convection and radiation terms. 

If Equations A-14 are used to reduce the equations to 
the form of Equation 2-9, the sensitivity equations become, 


[A'1{S3 4} + {dy} = 0 (B-9) 


where tha [A] matrix for the sensitivity equations is the 
same as the [A] matrix for the temperature equations, A-15. 


The {d} vectors for each parameter are listed as follows, 


n-1 
(beer Tart? + RCX, 8,4 7; 
At BBB, 


{d}1 = (rex, Sia" 
At /nee, 


(a-a1)h_o¢(T,~-U,) + RCX, S; 47-1 
f 1 = 1 ’ 
( re aw re 1,2 /enp, 


(4}2 = (frex, s, ."? 
At : /spp, 
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“yt 
. 


aS 


= 7 ae 
(« a2)hi.¢(Ta,-U,) + RCX, S, 3% y 
At BBB, 


~u,4 


_RP_(U,-U,)+eo0(U,4 jhe he (Tc 0) 
1 (0g 0,) +e0( 0, Ug) Mparbret Taw V1) 
( c Cc c 


~-1 
+ RCX,S; 4” 
At BBB, 


é. 6, At BBB 


RP, (U,-U,)+RCX;S; 57> 
by At sep, 


{d}5 = RM, U;_4-(RM,+RP,)0,+RP{U3,,+RCX;S; 5° 
6, 6, by At’ BB; 


(B-10) 
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APPENDIX C 


The HEATEST program follows. 
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